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ABSTRACT 

Magain et al. (2005) argued that the host galaxy of the quasar in HE0450-2958 is 
substantially under-luminous given the likely mass of its nuclear black hole. Using 
kinematical information from the spectra of the quasar and the companion galaxy, an 
ultra-luminous infrared galaxy, we test the hypothesis that the black hole powering 
the quasar was ejected from the companion galaxy during a merger. We find that the 
ejection model can be securely ruled out, since the kick velocity required to remove the 
black hole from the galaxy is ^ 500 km/s, inconsistent with the presence of narrow 
emission line gas at the same redshift as the quasar nucleus. We also show that the 
quasar in HE0450-2958 has the spectral characteristics of a narrow-line Seyfert 1 
galaxy and calculate a mass for its black hole that is roughly an order of magnitude 
smaller than estimated by Magain et al. The predicted luminosity of the host galaxy 
is then consistent with the upper limits inferred by those authors. 



1 INTRODUCTION 



HE0450-2958 is a bright quasar at redshift z = 0.285. 
HST images revealed that the system is double, with an 
ultra-luminous infrared galaxy (ULIRG) situat ed ~ 1.5 arc- 
sec fr om the qu asar | Bo yce et al 1 996: Ca nalizo fc Stockton! 
I2OOII) . Recentlv. lMaeain et alTBoOSh reported that the host 
galaxy of the quasar is substantially under-luminous, based 
on the quasar's luminosity and on a likely value for M., the 
mass of its nuclear supermassive black hole (SBH). Magain 
et al. proposed either that the quasar host galaxy is dark, 
or that an otherwise "naked" SBH had acquired gas while 
moving through intergalactic space. 

Here we examine these hypotheses in light of additional 
evidence from the spectra of the quasar and the compan- 
ion galaxy. The quasar spectr um reveals it to be a typ i- 
cal narrow- line Seyfert 1 galaxy iOsterbrock fc P ogge 1985), 
not a giant elliptical galaxy as assumed bv iMagain et aL 
1(200^. We infer a much smaller luminosity for the host 
galaxy, consistent with the upper limits derived by those 
authors. We also critically examine the most natural model 
for a "naked" SBH, namely, a SBH that was ejected from 
the com panion galaxy durin g the merger that created the 
ULIRG (iMerritt et al.ll2004) . We show that the ejection 
model can be securely ruled out, since the quasar spectrum 
indicates the presence of narrow emission line gas extending 
out to a distance of ~ 1 kpc from the nucleus that is moving 
at the same velocity as the broad-line gas. The narrow-line 
gas could not have remained bound to the SBH if it were 
ejected from the companion galaxy. 



2 SPECTRAL ANALYSIS 

HE0450-2958 was observed during November 27 2001 us- 
ing the UV Focal Reducer and low- dispersion spectrograph 
(FORSl) on Unit Telescope 1 of the VLT (PI: M. Courbin). 
The instrument was operated in MOS mode with a long-slit 
position angle of ~ 55° . This allowed the contributions from 
the ULIRG, quasar and the nearby G-type star to be gath- 
ered simultaneously in slit no. 9. In total, five spectra were 
obtained of HE0450-2958: three of 1200 s duration using the 
600B grism centered at 4620A, and two of 1800 s duration 
using the 600R (6270A) and 6001 (7940A) grisms. These 
data, including all relevant calibration files, were retrieved 
from the VLT data archive ^ . 

The data were reduced using standard IRAF^ routines. 
Bias and flat-field subtraction were carried out before wave- 
length calibration with HeArNe arcs. Cosmic ray subtraction 
was facilitated with a median combine, in the case of the 
three 600B exposures, and with a median fllter and rigorous 
visual inspection in the 6QQR and 6QQI cases. Background 
subtraction was performed by removing a third-order poly- 
nomial fltted to the sky components in the spatial direction. 
Flux calibration was carried out by fltting a 5700K (G-type) 
black body spectrum to the observed star, and scaling to 
the fluxes observed through the High Resolution Channel of 



^ "http : //archive . eso . org/" 

IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Radial velocities obtained from the peak wavelengths 
of the strongest emission lines along the slit. The locations of the 
galaxy, quasar and "blob" identified by Magain et al. (2005) are 
indicated. 



the Advanced Camera for Surveys aboard the Hubble Space 
Telescope (no. 10238, PL Courbin). 

Figure shows the velocities derived from the peak 
wavelengths of the strongest emission lines along the slit; 
the regions dominat ed by the ULIRG, the quasar and the 
iMaeain et alj (|200^ "blob" are indicated. In the region 
dominated by the galaxy, the velocities show the pattern of a 
rotation curve, with peak-to-peasuch a k amplitude of ~ 140 
km s~^. This value seems small for such a luminous galaxy, 
but the HST image suggests that the orientation of the slit, 
aligned to include the star, quasar and ULIRG, is far from 
the major axis of the galaxy. According to the rotation curve, 
the systemic velocity of the galaxy is 73950±20kms~^. In 
the region dominated by the quasar, the velocities obtained 
from the Ha, H/3, [NII]A6584 and [SII]A6717 emission lines 
do not vary within the uncertainties, and the average value 
is 73920±20 kms^^, indicating a blueshift relative to the 
systemic velocity of the galaxy of only 30kms~^, consistent 
within the errors with zero. The [OIII]A5007 emission line is 
blue-shifted relative to the other emission lines by approxi- 
mately 60 km s~^ . Howe ver such a blueshift in [OIII1A5007 is 
often observed in AG N llNelson fc Whittldll995l:lBian et all 
l2005l:lBorosor3l2005h . In the region dominated by the emis- 
sion of the blob, this blueshift is not observed, and the ve- 
locities derived from the [OIII] emission line are the same as 
those obtained from the other lines. In summary, our mea- 
surements show similar systemic velocities for the ULIRG, 
quasar and blob. 

The blue spectrum of the galaxy (Figure I^J shows a 
number of absorption-line features, in particular, high-order 
Balmer lines indicative of an intermediate-age (~ 10*yr) 
stellar population. In order to better quantify the age of the 
stellar population and to obtain an estimate for the stellar 
velocity dispersion, we performed a spectral synthesis using 
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Figure 2. The blue spectrum of the ULIRG, together with the 
synthetic spectrum and the residual between the two. The syn- 
thetic spectrum is ploted twice: superimposed on the galaxy spec- 
trum and shifted vertically for clarity. 



the code STARLIGHT of ICid Fernandes et all J2005l) . This 
code uses a basis of stellar population templates, each one 
corresponding to a given metallicity and age, to synthesize 
the galaxy spectrum, and it gives as output the contribution 
of each template to the total light at A4020A; the internal 
reddening; and the velocity dispersion. The contribution of 
the quasar, which contaminates both the continuum and 
the emission lines of the galaxy spectrum, has also been 
included in the synthesis. The results of the synthesis (Fig- 
ure 2) yield the following contributions to the light of the 
ULIRG at A4020: ~25% from the quasar continuum, ~50% 
from stars with ages around 10^ yr and ~25% from older 
stars. The total stellar mass is ~ 6 x 10^" M©, with the ma- 
jor uncertainty due to the quasar continuum contribution. 
The internal reddening is Ay—OA mag and the velocity dis- 
persion is (T = 190 ± 25 kms^^. The synthesis strongly 
points to a major burst of star formation ~ 10* yr ago. 



3 EJECTION HYPOTHESIS 

The low luminosity of the quasar host, coupled with its prox- 
imity to a ULIRG, leads naturally to the hypothesis that 
the SBH powering the quasar was ejected from the ULIRG 
following a merger. Two ejection mechanisms have been dis- 
cussed: gravitati onal radiation reco il during the coalescence 
of a binary SBH llFavata et all2004r) : or a gravitational shng- 
shot involving three SBHs, if the merger happened to bring 
a third SBH in to the center of a ga laxy containing an unco- 
alesced binary jMikkola fc Valtonen.l99Ql 

The quasar is displaced 1.5 arcsec from the center of the 
ULIRG, corresponding to a projected separation of ~ 6.5 
kpc. This is much greater than the galaxy's half-light ra- 
dius implying an ejection velocity Vkick comparable to the 
central escape velocity from the galaxy Vesc- The distri- 
bution of mass around the companion ga l axy o f HE0450- 
2958 is unknown. However, iTacconi et al.l 12002*) find that 
the light distributions in a sample of 18 ULIRGs are rea- 
sonably well fit by de Vaucouleurs profiles with effective 
(projected half-light) radii of _Re ~ 1 kpc, similar to those 
of luminous E galaxies. They derive kinematical masses 
of 0.3 X 10" Mr < M < 5 X 10"M(7,, consistent with 
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the estimate presented above from population synthesis of 
0.5-0.8 X 1O"M0. The mean stellar velocity dispersion in 
their sample is 180 km s~^, consistent with the population 
synthesis estimates presented above. We accordingly mod- 
elled the baryonic mass distribution around the ULIRG as if 
it were a normal, spherical E galaxy of mass Mgai, and used 
the empirical correlations between E -galaxy mass, luminos- 
ity, effective radius and Sersic index iMagorrian etaPllQQSl : 
iGraham fc Guzmanl20o3lFerrarese et al.l2005ll to derive its 
gravitational potential. 

Table Ogives the parameters of two mass models for 
the companion galaxy. Model 1 has a baryonic mass of 
1.5 X 10^°Mo, 

roughly the mass of a Mb « — 18 dE galaxy 
and a factor of ~two sma^er than the smallest ULIRG mass 
inferred by Tacconi et al. Model 2 has Mgai = 1.5 x lO" A^q, 
the stellar mass of a Mb ~ — 20 E galaxy, and close to the 
average mass of the galaxies in the Tacconi et al. sample. 
We also included a dark-matter halo; because the contribu- 
tion of dark matter to the gravitational force on scales ^ 10 
kpc is probably much less than that of the baryons, we con- 
sidered only a single halo model. As templates for the dark 
matter, we considered the four "galaxy-sized" halos in the 
iDiemand et alj ||2004) ACDM simulations, which have virial 
masses in the range 1.0 x 10^^ Mq ^ Mpm 2.2 x lO^^Mp . 
Fits to p(r) for these halos are given in lGraham et alj l)200!Tfl : 
based on these fits, central escape velocities lie in the range 
480 km s"^ < Vesc ^ 600 km s"\ and the in climb- 
ing to 10 kpc is 210 km s"^ ^ ^ 310 km s"\ Our 
adopted halo model (Table had a mass of 2.0 X IO^^Mq 
(virial mass 1.1 x lO^'^M©), half-mass radius 200 kpc and 
ce ntral escape velocity ~ 500 km s~^, similar to model Gl 
in iDiemand et al.l (|2004|) . By comparison, the virial mass 
of the Milky Way halo is believed to be 1 - 2 x lO^M© 
ijKlvDin et alJl2002h . 

In de Vacouleurs-hke mass models, the AF in climbing 
from the very center out to a distance of a few pars ecs can 
be con siderable due to the high nuclear density (e.g. lYound 
il97(fl l. The mass distribution of the companion galaxy fol- 
lowing the merger is unknown on these small radial scales, 
and in any case, the SBH would carry with it the mass of 
the inner few parsecs, modifying the potential. Accordingly, 
we placed the SBH initially at a distance of 10 pc from the 
galaxy center when computing post-ejection trajectories. 

Figure 121 shows the results. As expected, kicks of ~ 500 
km s~^ (Model 1) and ~ 1500 km s~^ (Model 2) are required 
in order for the SBH to climb a distance of 10 kpc from 
the center of the galaxy. This result essentially rules out 
radiation recoil as the origin of the kick, since the maximum 
ampl itude of the recoil i s believed to be less than 500 km 
s"^ llMerritt et alJl2004) and probably no more than 250 
km s~^ (lBlanjhe^^iri20o3) . Three-body recoils might still 
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Figure 3. Trajectory of a kicked SBH in two models for the mass 
distribution of the companion galaxy (Table 0. Black (thick) 
curves: Model 1; blue (thin) curves: Model 2. Curves are labelled 
by Vkick in km s'^. 



work, however the largest AV in a three-body interaction is 
experienced by the smallest body. But Figure constrains 
any ejection model in a number of other ways: 

1. A large Vkick implies a large velocity, V <; 300 km 
s~^, as the SBH moves past its current position. This is 
hard to reconcile with the essentially zero radial velocity 
difference between quasar and galaxy, unless the ejection 
velocity is fine-tuned or nearly perpendicular to the line of 
sight. 

2. The time for the kicked SBH to reach a distance of 
10 kpc is much shorter than 10** yr, again unless the gravi- 
tational potential and kick velocity are finely tuned. But the 
starburst occurred ~ 10* yr ago; thus, either the true sep- 
aration of the SBH from the ULIRG is much greater than 
10 kpc, or the ejection was delayed until a time of ~ 10* yr 
after the starburst. 

3. A large Vkick implies that the ejected SBH will carry 
very little mass with it as it departs the galaxy. Material 
orbiting the SBH with velocity v ^ Vkick before the kick 
will experience the kick as an adiabatic perturbation and 
will "instantaneously" acquire the specific momentum of the 
SBH. This argument suggests that the SBH will carry with 
it the mass contained initially within a region whose size 
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is less than r^ff, the radius at which the orbital velocity 
around the SBH is equal to Vkick, or 



GM. / Vkick \ 



-2 I'Vi 

lOMgo-aoo 



^ Vkick ^ 



pc 



(1) 



with Mg = M./10**Mq and 0-200 = O-/200 km s"\ Since 
Mg « 1 (see below) and Vkick 3cr (Figure I^J, the en- 
trained region will be of order 1 pc or less in size. This is 
probably sufficient to include the broad-line region (BLR) 
gas, which is expected to have a size ~ 0.3 pc based on 
the empirical scaling rel ation between BLR size and 5007A 
luminosity ^as m ct aObOOSl : iGreene fc Hcll2005^ ■ but not 
larger structures. 



4 SIZE OF THE NARROW-LINE REGION 

The radius of the narrow-line region can be estimated from 
the ionization parameter and density of the emitting gas and 
the ionizing luminosity of the QSO. For an ionizing photon 
luminosity, Q, the ionization parameter in gas of density 
n at a distance r from the source can be defined as [/ = 
Q/{'iTvR%i^iinc). We estimate the ionizing photon luminosity 
by extrapolating the slope o f the far UV (~ 6 00 - 1200A) 
continuum as determined bv lScott et alJ (|20o3) from FUSE 
observations. After corrections for interstellar extinction and 
absorption, Scott et al.'s power-law (/^ cx z^~") fit yields a 
spectral index a = —1.2 ±0.1 and a flux at lOOOA of 5.25 x 
IQ-^'^ergs"^ cm"^ Hz'^ HE0450-2958 has a steep soft X- 
ray continuum, wit h a photon index F = 3.1 in the 0.1-2 
keV ROSAT band llBrinkmann et al.lll997[) . Therefore, we 
extrapolate the UV power-law to a high energy cut-off of 0.1 
keV. Adopting a luminosity distance of 1458 Mpc (assuming 
Ho — 71 kms~^Mpc~^ and a matter density parameter 



flm = 0.27) the integration yields Q ^ 1.2 x 10'' photons 
s-\ 

The gas density was obtained from the relative inten- 
sities of the [SII]AA6717, 6731 doublet. Although partially 
blended in the spectrum, Gaussian fits to the lines are well 
constrained and yield a ratio leyn/Iersi — 1.01 ±0.05. This 
corresponds to a density n ~ 1000 cm~'^. 

We determined the ionization paramete r from the ra- 
tio o f the [OII]A3727 and [OIII]A5007 lines jBaldwin et alJ 
Il98ll) . The intensity ratio measured from the spectrum is 
^3727/^5007 ~ 0.13 ± 0.05. Using the empirical relation be- 
tween this ra ti o and the ionization parameter given by 
IPenston et"ai] (Il99dl . we obtain U ^ 0.014. This value 
is broadly consistent with the mo r e rece nt photoionization 
models presented bv lGroves ilaOiooi). 

Combining estimates for Q, n and U, we arrive at a ra- 
dius Rnlr ~ 1.5 kpc. This of course represents an ill-defined 
spatial average, since we have used integrated fiuxes of emis- 
sion lines representing different ionization states. Moreover, 
the determination of U is model-dependent and extrapo- 
lating the far UV power-law is, at best, a crude representa- 
tion of the EUV continuum.^ Based on these considerations, 
Rnlr may be uncertain by a factor ~ 2. It is, nevertheless, 
unsurprising that the NLR in a relatively luminous quasar 

^ However, the Ha photon luminosity is a: 10^® photons s~^ 
which, assuming photoionization, implies that our value for Q is 
underestimated by at least a factor ~ 3. 



extends to a few kpc. Direct imaging of bright, low-redshift 
quasars in [OIII] w ith HST reveals NLR sizes ranging from 
1.5 kpc to 10 kpc teen nert et"al]l2002h . Indeed, the [OIII] 
luminosity of HE0450-2958 (L[oj/j] « 3.6 x lO^^ergs"^) 
makes it comparable with the most luminous object in Ben- 
nert et al.'s sample, for which they determine an NLR radius 
of 10.5 kpc. 

Such a large size for the NLR rules out the possibil- 
ity that the NLR gas would remain bound to the SBH after 
ejection from the ULIRG (cf. equation 0. Post-ejection ac- 
cretion of the NLR gas from a cloud is also unlikely, since 
the radius of the Bondi accretion column, Vacc, is given by an 
equation similar to eq.0 after replacing Vkick by the relative 
velocity between SBH and gas cloud, implying Vacc <C 1 kpc 
unless the ejected SBH has fortuitously matched velocities 
with the cloud. 



5 MASS OF THE BLACK HOLE AND 

IMPLICATIONS FOR THE HOST GALAXY 
LUMINOSITY 

As shown in Figure 21 HE0450-2958 exhibits characteris- 
tics which unambi g uously identify it as a narrow-line Seyfert 
1 (NLSl; iGrupd (120041 ^'). Specifically, its broad Balmer 
lines have FWHM's ~ 1300 km s~^ (the conventional defi- 
nition requires FWHM < 2000 kms^^; JDstcrbrock & Pok3 
1^8^), it has strong optical Fell emission features and, 
as already noted, it has a steep soft X-ray photon con- 
tinuum. The currently-favored picture of NLSl's is that 
they represent an extreme AGN population character- 
ized by relatively low-mass SBHs bu t high accretion rates 
JPeterson et al.M2000t | BorosonI I 



tially super-Eddingtori ilBoU^ 



2002 ) — possibly substan- 



2003). It follows that esti- 



mating M, from the quasar luminosity ass uming a sub- 
Eddington accretion rate, as was done by iMaeain et alJ 
12005)), is likely to be misleading. Here we adopt what we 
consider to be a more robust approach, and estimate a virial 
mass based on the velocity dispersion {v) and radius (Rblr) 
of the broad-line region : Mbh ~ v^Rblr/ G iWandel et alJ 
ll999l:lKaspi et alJl2000l:IVestergaarJ2002^ . 

In this method, the BLR velocity dispersion is derived 
from the broad-line widths while the BLR radius is inferred 
from the radius-luminosity relation derived from reverbera- 
tion mapping (Kaspi et al. 2000; Peterson ct al. 2004). Here, 
we use the recent re vision of [K aspi ct al.| (|2000|)'s virial 
formula presented by iGreene Hq L200S] . This requires 
measurements of the H/3 FWHM {v = a/(3)/2x FWHM) 
and the continuum luminosity at 5100A(AL5ioo). Our mea- 
surements of these quantities yield FWHM(H/3) « 1270 
kms~^ and ALsioo ~ 4.6 x lO^^ergs"^, respectively. In- 
serting these values into equation (5) of Greene & Ho, we 
obtain Mbh = (9 ± 1) x 10^ AIq. Greene & Ho's alternative 
virial formula, which employs the luminosity and FWHM 
of the broad Ha line yields a consistent result, albeit with 
greater uncertainty: AIbh = (6!l^3) x 10^ Mq. These masses 
are subject to a systematic uncertainty of a factor ~ 3 re- 
lated to the poor ly-known structure, kinematics and aspect 
of the BLR fe.g.. lOnken et alJ i2004i ). Nevertheless, at face 
value, the virial method yields a SBH mass that is an order 
of magnitude less than the value M, « 8 x 10* Mq adopted 
bv lMaeain et alJ ^M^. 
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Figure 4. Comparing the quasar spectrum with that of 
PG1211+143, a radio-quiet quasar with steep X-ray spec- 
tral index that is classifi ed spectroscopically as a NLSl 
iConstantin fc Shieldj l2003l) . HE0450-2958 is plotted in ab- 
solute flux units, PGf2ff+143 in absolute units minus 
10~^^ergs~^A~^cm~^. Residuals are plotted along the bottom 
with a solid straight line highlighting zero. 

The host galaxies of NLSl's are spirals, often barred 
JCrenshaw et aLll2003l) . but relatively little is known about 
their systematic properties. The tight correlation between 
bulge velocity dispersion and SBH ma ss that characterizes 
quies cent elliptical galaxies and bulges (iFerrarese fc Merritd 
1200(3) also appears to b e valid for the bulges of active 
galax ies, including NLSls (|Ferrarese et all200lHBotte et all 
lioO^. Adopting M. = 9 x IO'^Mq, we infer a bulge veloc- 
ity dispersion cr « 180 km s~^. Near-IR bulge luminos ities 
also correlate tightly with M. l)Marconi fc Hun'3l200,^) : we 
infer a Jf-band absolute magnitude for the stars in the bulge 
of Mk ~ —23.4. Visual bulge ma gnitudes are more poorly 
correlated with M,. Adopting the IFerrarese fc FordI ()20o'il 'l 
relation gives an absolute blue magnitude Mb = — 18.9±0.5; 
alter nately, applying a B — K color correction of 4.0 to 
Mk llPeletier fc Balcellsl Il99d) gives Mb ^ -19.4 Com- 
puting Mb directly from a via the iFaber-JacksonI (119761) 
relation gives a similar value. Even more uncertain is the 
predi cted total (bu lge-fdisk) luminosity. Assuming an Sa 
host iW hittld 11 9921) implies a d isk-to-bulge ratio of ~ 1.5 
iSimien fc de Vaucouleur3ll986l) and a total visual magni- 
tude Mv ~ —21. While very u ncertain, this estimate is 
2.0 — 2.5 magnitudes fainter than lMaeain et all (120051) 's es- 
timate (—23.5 ^ Mv ^ —23.0) based on a ~ lOx larger 
assumed value of M,, and consistent with their conclusion 
that the host galaxy must be at least 4—5 magnitudes fainter 
than the quasar {My = —25.8). 



6 CONCLUSIONS 

The HE0450-2958 system consists of a ULIRG that expe- 
rienced a major starburst ^ 10* yr ago, situated at ~ 7 
kpc projected separation from a quasar having the spectral 
characteristics of a narrow-line Seyfert 1. The quasar host, 
presumably an early-type spiral galaxy, is expected to have 
a bulge luminosity Mk ~ —23.4 (M y ~ —21). Upper li mits 
on the luminosity of the quasar host (jMaeain et alJ2005l) are 
consistent with the expected total luminosity for a galaxy 



of this type. The SBH that powers the active nucleus ap- 
pears to be accreting at a super-Eddington rate, L/Le ~ 3, 
similar to the accretion rates inferred in other NLSls. Ejec- 
tion of the SBH from the ULIRG is unlikely for a number 
of reasons, the strongest of which is the presence of narrow 
emission line gas at the same redshift as the quasar nucleus; 
this gas could not have been retained if the SBH was ejected 
from the companion galaxy. We find no compelling evidence 
that the quasar in HE0450-2958 is either a "naked" SBH 
ejected from its host galaxy, or that it has an anomalously 
dark host galaxy. 
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